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ABSTRACT 

X-ray emission from about ten protostellar jets has been discovered and it appears as a feature 
common to the most energetic jets. Although X-ray emission seems to originate from shocks internal 
to jets, the mechanism forming these shocks remains controversial. One of the best studied X-ray jet 
is HH 154 that has been observed by Chandra over a time base of about 10 years. We analyze the 
Chandra observations of HH 154 by investigating the evolution of its X-ray source. We show that 
the X-ray emission consists of a bright stationary component and a faint elongated component. We 
interpret the observations by developing a hydrodynamic model describing a protostellar jet originating 
from a nozzle and compare the X-ray emission synthesized from the model with the X-ray observations. 
The model takes into account the thermal conduction and radiative losses and shows that the jet/nozzle 
leads to the formation of a diamond shock at the nozzle exit. The shock is stationary over the period 
covered by our simulations and generates an X-ray source with luminosity and spectral characteristics 
in excellent agreement with the observations. We conclude that the X-ray emission from HH 154 is 
consistent with a diamond shock originating from a nozzle through which the jet is launched into 
the ambient medium. We suggest that the physical origin of the nozzle could be related to the 
dense gas in which the HH 154 driving source is embedded and/or to the magnetic field at the jet 
launching/collimation region. 

Subject headings: Hydrodynamics; Herbig-Haro objects; ISM: jets and outfiows; X-rays: ISM 



1. INTRODUCTION 

Thanks to the capabilities of current X-ray observa- 
tories (Chandra and XMM-Newton), a new class of X- 
ray sources has been discovered in the past few years, 
namely those associated to protostellar jets. This class 
consists of about ten members, HH 154 being one of 
the best studied. The X-ray source associated with 
HH 154 has been p reviously observed by XMM-Newton 
iFavata et al."'2002) and by Chandra (iBallv et al.ll2QQ3l : 
iFavata et al. 2006). These observations had shown that 
the X-ray source is stationary over a period of 4 years 
with a possible transie nt and faint compon ent in the 2005 
Chandra observations (|Favata et al.ll2006[ ). 

Hydrodynamic models of the interaction between a 
continuously ejected supersonic jet and the ambient 
medium predict X-ray emission caused by shocks at the 
interaction front (|Bonito et al.|[2QQ3 . 120071 ). X-ray emis- 
sion close t o the base of the jet i s predicted in the case of a 
pulsed-jet CBonit o et aTll2010bl laf). The pulsed-jet model 
is very effective in reproducing the features and variabil- 
ity observed in most of the X-ray emitting jets from low- 
Electronic address: sbonito@astropa.unipa.it 

^ INAF - Osservatorio Astronomico di Palermo, P.zza del Par- 
lamento 1, 90134 Palermo, Italy 



mass young stellar objects ("Bonito et al.J l20l"Qal ). In the 
pulsed-jet scenario, iBonito et al. (20108]) provided pre- 
dictions on future observations of HH 154 and claimed 
that a stationary X-ray source on a period of « 10 years 
is unlikely. For a stationary X-ray source, they suggested 
that a nozzle formed by either a dense medium or a mag- 
netic field may be at work at the launching site, leading 
to the formation of a stationary diamond shock at the 
base of the jet. 

The new Chandra observation of HH 154 (collected 
in 2009) provides a time base of 8 years to analyze 
the spectral and morphological variability of the X-ray 
source. Here we analyze the 2009 observations and com- 
pare the results of the data analysis with the p revi- 
ous Chandra (Ballv et al. 2003; Fa vata et al.ll2006h and 
XMM/Newton (;Favata et al.i .2002i) observations. We 
have developed a model describing a protostellar jet out- 
fiowing from a nozzle and compare the X-ray emission 
synthesized from the model with the Chandra observa- 
tions. Finally we interpret the new and the previous ob- 
servations in the light of our model results. In Sect. 2 we 
describe the Chandra observations and their analysis; in 
Sect. 3 we discuss the hydrodynamic model of jet/nozzle 
and the comparison between model predictions and ob- 
servations; in Sect. 4 we discuss the results and draw our 
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conclusions. 

2. CHANDRA OBSERVATIONS OF HH 154 

We present the new Chandra/ACIS-S data of HH 154 
(PI Schneider; ObsID 11016; texp = 65.2 ks) centered 
at (04:31:34.998, +18:07:51.95) (FK5), performed on 29 
Dec. 2009, th at provides 8 years of time base between 
the first f2001. iFavata et al.ll2002l ) and the last Chandra 
observations of HH 154, w i th th e intermediate observa- 
tion of 2005 (Favata et al.l 120061 ). To derive a uniform 
comparison between the three Chandra data-sets, we re- 
processed all the data with the same method, using the 
latest CI AO 4.3 package, and filtered the data to restrict 
the energy to the 0.3 — 4.0 keV band, as in Fav ata et al.l 
(^2006). Events were extracted for all observations from 
regions around the source and the background, near the 
position of the HH 154 je t (4:31:34.10, +18:08:04.9), fol- 
lowing B allv et all (|2003h . We defined a box 3.5'' x 5'' 
in the 2001 data-set, a box 3.5'' x 5.5" in the 2005 data- 
set, and a box 3" x 6" in the 2009 data-set^, to extract 
the source events, and four 10 pixel radius circles for the 
background, from source-free parts of the image around 
the source. We extracted the spectra of the three data- 
sets and produced the ancillary response files, arf and 
rmf, by using SPECEXTRACT. We grouped the spec- 
tra to have a minimum of 5 counts per bin. 

We also applied the sub-pixel repositioning algorithm 
available in CIAO 4.3 (EDS ER) to the Cha ndra images 
to refine the event positions (jLi et al.ll2004f ). We resam- 
ple the improved images at 0.25" pixel size, obtaining 
images with one-half of the native ACIS pixel scale. Fig. 
[1] shows the X-ray source associated with HH 154 in 2001 
(first panels), 2005 (second panels), and 2009 (third pan- 
els), and the X-ray maps synthesized from our hydrody- 
namic model (discussed in Sect. [3]) with the same spa- 
tial resolution as ACIS (last panels). The upper panels 
show the images with the native ACIS pixel size (0.5"); 
the central panels show the images with sub-pixel repo- 
sitioning algorithm applied and resampled with a pixel 
size 0.25"; the lower panels show the improved resampled 
images smoothed using a Gaussian kernel of width 0.5". 

2.1. Results 

The three Chandra data-sets show that the X-ray emis- 
sion is mainly located at the base of the HH 154 jet in 
all epochs, near the driving source. iBallv et al.l (|2003l ) 
found that the X-ray source is displaced by 0.5" — 1" 
with respect to the L1551 IRS5 driving source. The X- 
ray source consists of a bright knot which appears to 
be stationary in the time period covered by the obser- 
vations (in the following the "stationary" component) 
and an elongated structure directed away from the driv- 
ing source (in the following the "elongated" component) 
showing variability in the three data-sets (see Fig[T]). In 
particular, the latter component appears as a knot in 
images with EDSER applied (central and lower panels 
in Fig. [T]) which is much fainter than the knot of the 
stationary component. The position of the faint knot 
appears to be different in the three data-sets. In partic- 
ul ar, for the 2001 and 2005 data, we confirm the results 
of iFavata et al.l (|2006li ). i.e. an elongation of the X-ray 

^ The size of the boxes is different for different years to account 
for the variation of the morphology of the X-ray source. 



source corresponding to a proper motion of 0.7"/yr ap- 
proximately westward away from the driving source of 
the HH 154 jet. As for the 2009 observations, there is 
an hint of a faint knot closer to the brighter station- 
ary source than in 2005. Note that recently an asym- 
metry in the Chandra PSF has been discovered^, lo- 
cated at P.A. = 195 — rollangle{^25y ^ corresponding 
to P.A. = 269.4° in 2009 observations, i.e. roughly the 
direction of the extension of the X-ray source detected in 
HH 154. We checked if this instrumental effect may influ- 
ence the observed morphology of the X-ray source and 
found that the asymmetry of the PSF does not affect 
our images on scales larger than 1 arcsec. The elongated 
structure visible in the images, therefore, is not an ar- 
tifact of the instrument. The maximum length of the 
whole X-ray source is ^ 5" (^ 700 AU at i:> = 140 pc). 

We performed a spectral analysis of the individual data 
sets^. All spectra are well fitted by an absorbed thermal 
plasma (APEC in XSPEC). The count rates of the three 
observations (0.76 ± 0.10 cnts/ks in 2001, 0.65 ± 0.08 
cnts/ks in 2005, and 0.89 ± 0.12 cnts/ks in 2009) are 
compatible within the Poisson error, indicating that the 
X-ray luminosity is constant. Moreover the source shows 
no spectral variability in terms of temperature (T) and 
emission measure {EM) in the timescale analyzed here. 
We have thus fit the three data-sets simultaneously find- 
ing the values reported in Table [H 

The Nil is well constrai ned (more than in the analy- 
sis of XMM-Newton data; iFavata et al"]l2002l ) by fitting 
the three data-sets simultaneously but it cannot be con- 
strained by the individual data-sets. In the latter case 
we fixed A^h to the value derived in the simultaneous fit- 
ting of the three data-sets. The three data-sets and the 
best-fit model are shown in Fig. [21 

To investigate possible spatial variations of the spec- 
tral properties, we selected two regions, confining the 
stationary component (blue box in Fig. [T]) and the elon- 
gated component (green box in Fig. [T]), and computed 
the median photon energy, MPE^ within each region in 
2001, 2005, and 2009. MPE has been proved to be a 
robust indicator of the spec tral pr operties of a source 
in the case of low statistics (Hong et al.jl200j). In each 
data-set the total number of counts considered is ^ 50, 
in the stationary component, and ~ 10, in the elongated 
component. The stationary component shows no tem- 
poral variability in the three epochs, its MPE varying 
in the range MPE^ = 1.35 - 1.42 keV. The median en- 
ergy of the elongated component MPEq is always lower 
than MPEs. By considering the three observations alto- 
gether MPEs ~ 1.4 keV, while the elongated component 
appears softer with MPEq ^1.0 keV. 

This result is supported by the Kolmogorov-Smirnov 
test we performed to check if the photon energy distribu- 
tions of the two regions differ significantly and vary with 
time. We found that the stationary and the elongated 
component are compatible with being constant over the 
time base analyzed. We also found that the two compo- 
nents are dissimilar at the 99.99% confidence level. In 
the light of this result, we have fitted simultaneously the 

^ See |http : //cxc . harvard . edu/ciao/caveats/psf .artifact . html] 

for detciTIsl 

We ha ve verified that ou r pr ocedure reproduces the results 
obtained bv IBallv etHl (|2003l ) and IFavata et al.l (|2006h . 
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Fig. 1. — Upper panels: X-ray count rate maps (0.3 — 4 keV) of the 2001, 2005, and 2009 data-sets (first to third panels) at the native 
ACIS pixel size (0.5^^). The last panel shows the synthetic X-ray map of the base of the jet as derived from the model at the same spatial 
resolution as Chandra/ACIS. Central panels: The same as upper panels, but resampled at 0.25^^ EDSER technique has been applied on 
these data as explained in the text. Lower panels: The same as central panels, but with a smoothing applied on the images (Gaussian 
kernel of width 0.5^^). In each panel No rth is up and East is left. The angular size of each panel is ~ 7" x 7" . The P.A. of the PSF 
asymmetry (see explanation in Sect. 12. 1|) is 93° in 2001, 102° in 2005, and 269° in 2009. Note that this asymmetry can produce artificial 
extension on angular scales up to 1^^, as discussed in the text. 



TABLE 1 

Best-fit values derived from the three data-sets simultaneously fitted. 
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Fig. 2. — Best-fit X-ray spectral model superimposed on the 
three Chandra data-sets (2001 in black, 2005 in red, and 2009 in 
green) fitted simultaneously. 



data collected in the three epochs for each of the two 



components and found that the stationary component is 
described by a plasma with T = (7— 14) x 10^ K and the 
elongated component by a plasma with T < 7 x 10^ K 
(assuming the same Nh for the two spectra), confirming 
that the elongated component is softer. 

3. THE MODEL 

X-ray observations of HH 154 suggest the presence of 
a steady shock at the base of the jet over a time interval 
of about 8 years. Examples of quasi- stationary shocks 
in jets are the diamond shocks, namely stationary shock 
patterns appearing in supersonic jets when the jet ma- 
terial outflowing from a nozzle is slightly over or under- 
expanded or, in other words, when the pressure of the gas 
exiting the nozzle is below or above th e pressure of the 
ambient medium. iBonito et al.l (|2QlQaf ) suggested that a 
diamond shock forming near the launching/collimation 
site of the jet is the most likely mechanism leading to a 
stationary X-ray source in HH 154. 

To test the above idea, we extended the hydrodynamic 
model of iBonito et al.l (|2QQ7[ ) (to which the reader is re- 
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Fig. 3. — Enlargement of the nozzle through which the jet is 
ejected into the ambient medium. 

ferred for more details), describing the interaction be- 
tween a continuously ejected supersonic jet with the un- 
perturbed medium, by adding a diverging nozzle at the 
base of the jet included as an impenetrable body (see 
Fig. [3|). The model is described by the hydrodynamic 
equations solved in two dimensions (adopting cylindri- 
cal coordinates and assuming axisymmetr y with the jet 
axis coincident with the axis of symmetry; iBonito et al.l 
l2QQ7f ) and takes into account the thermal conduction (in- 
cluding heat flux saturation) and radiative losses from 
optically thin plasma. The calculations were performed 
using FLASH, a well tested a daptive mesh refinement 
multiphysics code (jFrvxell et a l. 2000). 

The jet with temperature Tj = 10^ K outfiows from 
the nozzle and propagates through an initial uniform 
ambient medium with the same temperature. The ini- 
tial jet radius depends on the nozzle size, chosen accord- 
in^ly with observations of the order of tens of AU (see 
iBallv et~al] [2003). We have explored a wide space of 
model parameters defined by the particle number den- 
sity of the jet ranging between nj = 3 x 10^ and 5 x 10^ 
cm~^, the density of the ambient medium ranging be- 
tween TT-a = 3 X 10^ and 10^ cm~^, the throat of the noz- 
zle with radius ranging between i?th = 15 and 200 AU, 
and the jet velocity ranging between = 1000 and 1500 
km s~^. We assume that the direction of propagation of 
the jet is perpendicular to the line of sight. Our best-fit 
model is characterized by nj = 3 x 10^ cm~^, na = 3 x 10^ 
cm-3, iij = 1500 km s-\ and Rth = 100 AU. The jet is 
overexpanded and travels through an initially denser am- 
bient medium (light jet scenario; see Bonito et al. 2007) 
with ambient- to-jet density contrast u = Ua/rij = 10. 

The computational domain is (1000 x 4000) AU 
7 X 27 arcsec dit D = 140 pc). Infiow boundary condi- 
tions are imposed at 2: = and r < i?th, axisymmetric 
boundary conditions are imposed along the jet axis (con- 
sistent with the adopted symmetry) , and outfiow bound- 
ary conditions are assumed elsewhere. The maximum 
spatial resolution achieved by our simulations is ~ 2 AU, 
using: 5 refineme nt levels with the PARAMESH library 
([MacNeice et al.l ^2000). which corresponds to covering 
the initial jet radius of our best-fit model with 50 grid 
points. The grid resolution is increased inside the noz- 
zle and around the nozzle exit to capture the diamond 
structure forming there. Note that the pixel size of the 
hydrodynamic model at its maximum spatial resolution 
(2 AU) corresponds to a spatial resolution of 0.014'^ at 
the distance of 140 pc, i.e. much lower than the nominal 
spatial resolution of ACIS (0.5''). As discussed below. 



this implies that the synthetic maps derived from the 
hydrodynamic model needs to be rebinned to match the 
ACIS pixel size in order to compare the model results 
with the observations. 

3.1. Synthesis of X-ray emission 

Following iBonito et al.l (|2007[ ) , we synthesized the X- 
ray emission associated with the jet from the model re- 
sults, recovering the 3-dimensional spatial distributions 
of mass density and temperature. We then derived the 
emission measure, EM^ and temperature, T, for each 
fiuid element and the distribution EM{T) integrated 
along the line of sight, for each element, in the tem- 
perature range 10^ — 10^ K (using 75 bins equispaced 
in logT). From the EM(T), we synthesize the maps 
of the X-ray emission and the focal plane spectra using 
the APEC spectral code, and considering photon count 
statistics comparable with that of the observations. The 
source is assumed to be at a distance D = 140 pc. We 
filtered the emission through the Chandra/ ACIS instru- 
mental response and the interstellar column density at 
the best-fit value Nh = l.2x 10^^ cm-^ (see Sect. Elj). 
In particular, to compare the synthesized images with 
the different observations, we have used the rmf and arf 
files generated for the three data-sets (2001, 2005, and 
2009) to account for the correction to the charge transfer 
inefficiency and the degradation in low-energy response 
due to contaminant buildup on the optical blocking filter. 

To directly compare the synthesized images with the 
observations, we rebinned the model images (whose spa- 
tial resolution is 2 AU, corresponding to 0.014'' at 140 
pc) so as to have the same bin size as the Chandra images 
as shown in Fig. [I](0.25", or 0.5"). 

We also convolved the synthesized X-ray image with 
the specific PSF, created at the proper energy for each 
data set, by using the Chandra standard analysis tools. 
We finally included Poisson fiuctuations to mimic the 
photon count statistics. 

3.2. Results 

Figure [H shows the spatial distribution of density in a 
slice in the (x, z) plane (left panel) for the best-fit model 
at time t ~ 120 yrs and a close up view of the nozzle 
site (right bottom panel). After the jet exits the nozzle, 
it propagates through the ambient medium forming a 
shock with temperature T ^ 10^ K at its head and a cold 
cocoon enveloping it. The cocoon is characterized by low 
temperatures (T < 10^ K), is dominated by the radiative 
cooling, and is strongly perturbed by the hydrodynamic 
instability developing there as the jet progresses through 
the medium, the thermal conduction being inefficient in 
damping the instability (see Bonito et al.ll2010al lbl). 

The nozzle determines a diamond-shaped shock past 
the nozzle exit with peak temperature T ^ 8 x 10^ K (see 
right bottom panel in Fig. H]). This diamond structure 
has its origin inside the nozzle and appears as a shock 
emerging from the nozzle and reflecting just past of the 
nozzle exit. After its formation (t ^ 50 yrs), the diamond 
shock is almost stationary until the end of the simula- 
tion for ^ 100 yrs. The thermal conduction is rather 
efficient in the post-shock region given the high temper- 
atures there (T > 10^ K) and is crucial in stabilizing the 
diamond structure, damping the hydrodynamic instabil- 
ity developing past the nozzle exit. Auxiliary simulations 
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performed without the thermal conduction have shown 
that the diamond shock would be unstable if the ther- 
mal conduction is neglected, the hydrodynamic instabil- 
ity heavily perturbing the flow structure at the nozzle 
exit. 

Analyzing the X-ray emission synthesized from the hy- 
drodynamic model, as described in Sect. 13. ![ we inves- 
tigated both the morphology and spectral properties of 
the synthetic X-ray sources. The X-ray emission from 
the modeled jet consists of two main features: a quasi- 
stationary source associated with the diamond shock at 
the jet base and a moving source associated with the 
shock at the head of the jet. The latter is a transient 
feature we are not interested in^ and does not influence 
the evolution of the diamond shock at the base of the 
jet; therefore we will not discuss its properties in the 
following. 

The X-ray luminosity of the diamond shock is Lx ^ 
5 X 10^^ erg and is stationary over ^100 yrs. This value 
is similar to that observed for HH 154 that is almost 
stationary in about 8 years. 

By comparing the total flux derived from the model 
with the speciflc rmf and arf response of each data-set, we 
have verifled that the degraded QE of the instrument in 
the time baseline analyzed affects the synthesized count 
rate for less than 7%. This conflrms that the source flux 
can be assumed constant over 8 yrs, within the Poisson 
errors. 

The right upper panel in Fig. |4] shows the synthetic X- 
ray emission arising from the shock integrated along the 
line-of-sight. Most of the emission originates just behind 
the shock in a bright and compact knot with temperature 
T 8 X 10^ K. The knot is surrounded by a diffuse re- 
gion elongated along the jet axis, characterized by lower 
temperatures {T ^ 1 — 2 x 10^ K). Figure [5] shows the 
proflles of density and temperature along the jet axis in 
the region where the diamond shock forms. 

We found that the spectrum synthesized from the hy- 
drodynamic model, as explained in Sect. 13.11 can be flt- 
ted with one isothermal component which is compatible 
with that derived from the three data-sets of Chandra. 

We rescaled the synthetic X-ray image shown in Fig. |4] 
to the Chandra/ ACIS pixel size (last panels in Fig. [1]). 
The emission within the nozzle is assumed to be totally 
absorbed. We found that the spatial scales of the X-ray 
emitting diamond shock at the same spatial resolution of 
Chandra are consistent with the size of the HH 154 X-ray 
emitting source: a synthetic X-ray source of a few arcsec 
at the base of the jet consisting of a bright point-like com- 
ponent surrounded by a faint and elongated component 
along the jet axis. 

In Fig. [6] we compare the smoothed 2001 image with 
a bin size 0.25'^ (left panel) with the X-ray source de- 
rived from the model at its maximum spatial resolution, 
0.014^', (right panel). The 2001 image contour is super- 
imposed on the modeled source. 

4. DISCUSSION AND CONCLUSIONS 

^ The X-ray emission at the head of the jet forms because we 
simulate the early evolution of the jet when its head is close to 
the stellar driving source. We expect that the head of the jet 
progressively cools down because of radiative cooling as it goes 
away from the driving source, eventually emitting mainly in the 
optical band. 



The analysis of the observations of HH 154 in three dif- 
ferent epochs with Chandra reveals a faint an d elongated 
X-ray source displaced by 0.5 — 1 arcsec (|Ballv et al.l 
I2OO3I ) from the LI 551 IRS 5 protostar (the driving source 
of the jet) along the jet axis. The source appears to be 
quasi-stationary over a time base of ^ 8 yrs without ap- 
preciable proper motion and variability of X-ray lumi- 
nosity and of temperature. The morphological analysis 
shows that the X-ray source consists of a bright station- 
ary component with temperature T > 7 x 10^ K sur- 
rounded by an elongated cooler component extended in 
the direction away from the driving so urce, with tem- 
peratu res T < 7 X 10^ K. Very recent Iv [Schneider et al.l 
(|201ll ) analyzed the same data-sets finding similar ob- 
servational results in terms of X-ray luminosity, spectral 
parameters, and morphology, independently showing the 
robustness of the derived parameters that form the basis 
of our comparison with a simulation of the jet based on 
detailed hydrod ynamic models. 

As shown in Bonito et al. (2008), the X-ray source is 
not perfectly aligned with t he optical j et observed in HH 
154 (see Fig. 13 in Bonito efail 120081 ). In fact the HST 
images of .Fridlund et al.. (.2005f ) show that the optical jet 
from HH 154 is along P.A. ^ 254° (see also iPvo et al.l 
I2OO2D , while from t he X-ray data we derive a P.A. ^ 
270°. iBonito et al.l (jlOlOa) suggested that an ejection 
direction varying in time could explain the misalignment 
between the X-ray source and the optical jet. Since the 
jet driving source, LI 551 IRS5, is known to be a binary 
system (Bieging fc CohenI Il985t ). a jet precession could 
be induced due to the presence of the companion star. 

The absorption column density derived from the anal- 
ysis of the three data-sets is too low if compared with the 
150 ma g of absorpt i on of L155 1 IRS 5, confirming the re- 
sults of Ba llv et atl (|200 3>) and'Favata et alj (|2006D . This 
fact together with the evident displacement of 0''.5 — 1'' 
of the source from LI 551 IRS 5 and the lack of temporal 
variation in the X-ray fiux and spectral properties sug- 
gest that the X-ray emission detected in the three epochs 
unambiguously arises from the jet and cannot be of stel- 
lar origin. 

The observations suggest therefore that the X-ray 
emission of HH 154 or iginates in a s t anding shock located 
at the base of the jet. iBonito et al.l (|2010al ), by analyzing 
the X-ray emission arising from a pulsed jet model, have 
discussed the possibility to produce a standing shock at 
the base of the jet as a result of multiple self-interactions 
of plasma blobs ejected in different epochs by the driv- 
ing source and concluded that this mechanism is unlikely. 
In addition these authors suggested that the most likely 
mechanism leading to a standing shock over a period > 5 
yrs might be a diamond shock forming near the launch- 
ing/collimation site of the jet. 

This idea is expanded here by developing a model of 
jet outfiowing through a nozzle. We found that, in such a 
configuration, a standing diamond shock forms just past 
the nozzle exit at the base of the jet. The shock is stabi- 
lized under the action of the thermal conduction which 
damps the hydrodynamic instability developing within 
the cocoon and heavily perturbing the fiow structure. 
We found that the X-ray emission arising from the dia- 
mond shock has morphology and spectral characteristics 
in excellent agreement with those derived from the three 
data-sets of HH 154. The model also predicts that the 
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Fig. 4. — Density map (left panel) with the 2 x 10^ K contour superimposed, an enlargement of the base of the computational domain 
(lower panel on the right), and the X-ray map synthesized from the model (upper panel on the right) at the maximum spatial resolution 
achieved from the simulation. 100 AU correspond to at 140 pc. 
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Fig. 5. — Profiles of particle number density (solid line) and 
temperature (dashed line) along the jet axis in the region where 
the diamond shock forms. 

X-ray emitting plasma of the diamond shock cools down 
at larger distances from the driving source as inferred 
from the observations. In fact, from the Chandra data, 
we found that the MPE increases toward the base of the 
jet. While we cannot rule o ut that this result is a ssoci- 
ated with variations in Nh (jFridlund et al.l l2QQ5l found 
an increasing absorption towards the driving source along 
the jet axis), higher values of MPE can be indicative of 
higher temperatures. Such variations of plasma temper- 
ature would be naturally explained by our model. We 
conclude therefore that HH 154 offers the first evidence 



Fig. 6. Smoothed 2001 image of the HH 154 X-ray source (left 
panel) with a pixel size of 0.25^^ compared with the high resolution 
(0.014^^) image derived from the model (right panel). The contour 
of the observed source is superimposed on the images. Since we 
verified that the direction of the PSF asymmetry in 2001 is along 
the counter jet and since this asymmetry can produce artificial 
structures only out to 1" (four pixels at the 0.25^^ pixel size of 
the image on the left panel), the observed elongation is not an 
artifact due to the PSF asymmetry. The angular size of each panel 
is ~ X A" . In both panels. North is up and East is left. 

of a standing diamond shock at the base of the jet prob- 
ably near the jet launching/collimation region. 

We can infer a characteristic size, Ighi of the X-ray 
emitting source from the spectral analysis, using the 
value of the EMbest-fu^ and from the hydrodynamic 
model results, using the maximum particle density value 
derived from the diamond shock modeled. In particular, 
we find Ish > = 4 X 10^^ cm, where V is the volume 
derived from the EM = O.Sn'^V (following iFavata et al.l 
[20020 and Umax ^ 10^ cm~^ (see the peak of density 
in Fig. [5]) at 140 pc. This value is in good agreement 
with both the observed radius of HH 154, rj ~ 30 AU 
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(see discussio n on this paramete r in both models a nd ob- 
servations in lBonito et al.|[2QQ7l and feonito et al.i r2QQ8). 
and with the size of the diamond shock modeled. There- 
fore, although the spatial resolution of the ACIS data, 
improved by sub-pixel techniques, is more than an order 
of magnitude lower than that achieved by our numerical 
model, its diagnostic power allows us to infer detailed 
information on physical scales comparable to numerical 
simulations. 

The data, after this analysis, clearly show the pres- 
ence of an elongated structure on the right side of the 
main source. Concerning an eventual evolution of this 
elongated structure the most conservative interpretation 
could be that it is steady and we can hardly constrain 
its features, due to the limited photon statistics. How- 
ever, given the insight provided by the model and the 
evidence that stellar jets flows are inherently variable, an- 
other interpretation is that we are observing the diamond 
shock variability and possibly knots forma tion and mo- 
tion d ue to the changes of the jet flow (cf. iBonito et al.l 
l2QlQal for example of effects of variable jet flows). The 
sequence of the smoothed images of HH 154 with a spa- 
tial scale of 0.25'' shown in Fig. [T] suggests the presence 
of subsequent knots with a detectable proper motion. In 
particular, by comparing the 2001 and 200 5 data-sets, 
we confirm the results of Favat a et al.l (j2006l ) who found 
a detectable westward proper motion of the elongated 
component of the X-ray source, corresponding to ~ 500 
km/s; we find a hint of elongation again westward in the 
2009 observations away from the jet driving source, but 
closer to the stationary source than in the 2005 obser- 
vations. This evidence may suggest the presence of a 
newly formed knot propagating away from the diamond 
shock. However note that, in the 2009 observation, the 
PSF asymmetry is directed almost along the X-ray ex- 
tension axis, and therefore it could influence the X-ray 
source elongation up to an angular scale ^ V\ There- 
fore, both the evidence of a standing shock at the base of 
the HH 154 jet over a 8 yrs timebase and a moving knot 
in 2005 together with a hint of a new elongation in the 
2009 smoothed image, indicate the scenario of a nozzle, 
creating the standing sho ck, in t he presence of a pulsed 
jet, as described in Bo nito et al.l (pOlOa) . which may ac- 
count for the moving/elongated component. No matter 
how one interprets the observations, there is a clear need 
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for future observations of HH 154. 

The physical origin of the nozzle could be related to 
the dense gas in which the LI 551 IRS 5 protostar is em- 
bedded and/or the intense stellar magnetic field at the 
jet launching/collimation region. In the hypothesis that 
a magnetic nozzle causes the diamond shock observed in 
HH 154, the Chandra observations and the comparison 
with our model offer the possibility to constrain the mag- 
netic field strength near the jet launching/collimation re- 
gion. In fact, the model provides the total plasma pres- 
sure {psh + Psh'^sh) X-ray emitting diamond shock 
that reproduces the observations, where Psh, Psh and u^h 
are the pressure, mass density and velocity in the post- 
shock region close to the nozzle exit, respectively. Then, 
assuming the plasma f3 = {p^h + Psh'^sh)/(^^/^'^) ^ 
where B is the magnetic field strength, we derive B ^ 5 
mG in the magnetic nozzle at the base of the jet. In- 
ter estingly, thi s valu e is consistent with that inferred 
by iBallv et al.l (j2QQ3f ). namely B = 1 - 4 mG, in the 
context of shoc k s asso ciated to jet collimation, and by 
ISchneider et al.l (|2Qllf ). who find B ^ 6 mG, which is a 
reasonable v alue at the jet basis n ear the driving source, 
according to lHartigan et al.l (|2QQ7f ). We suggest therefore 
that the comparison between our model and the X-ray 
observations of HH 154 may allow us to probe the launch- 
ing/collimation region near the driving source, very dif- 
ficult to be directly observed in systems so obscured. 
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